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The magnetic moment anomaly — 2)/2 of the positive muon has been measured at the Brookhaven 

Alternating Gradient Synchrotron with an uncertainty of 0.7 ppm. The new result, based on data taken in 2000, 
agrees well with previous measurements. Standard Model evaluations currently differ from the experimental result 
by 1.6 to 3.0 standard deviations. 


1. Introduction 


2. Experimental Setup and Data Analysis 


According to Dirac’s theory, the gyromagnetic 
ratio g of a lepton is exactly 2. Deviations from 
this prediction are caused by radiative corrections 
to the lepton-photon vertex due to quantum field 
fluctuations. The anomaly of the electron is cur¬ 
rently known at the level of 4ppb [1] and well 
in agreement with the Standard Model. Since 
the contribution of heavy virtual particles to the 
anomaly a = {g — 2)12 is proportional to the 
square of the mass scale, the sensitivity of the 
muon is enhanced by a factor ~ 40000 

relative to the electron. At the present level of 
precision, the muon anomaly probes QED, 
weak and hadronic contributions. 


The general technique of the experiment at 
BNL is the same as that of the precursor experi¬ 
ment at GERN [2]. Polarised muons are stored 
in a highly uniform magnetic dipole field with 
electrostatic quadrupoles [3] providing vertical fo¬ 
cussing. The muon spin processes relative to the 
momentum vector with the angular frequency 
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provided that P ■ B = 0. The dependence of iVa 
on the electric field E (second term in Eq. (1)) is 
eliminated by storing the muons at the “magic” 
7 ^ = 29.3, corresponding to a momentum p^ = 
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3.094 GeV/c. In this case, is given by simulta¬ 
neous measurement of to a and {B), the magnetic 
field averaged over the spatial muon distribution 
in the storage region. As explained later, Wa is 
reflected by the rate of decay positrons above a 
certain energy threshold. The magnetic field B is 
measured in terms of the free proton precession 
frequency Wp in this field using nuclear magnetic 
resonance (NMR) probes. Then, can be ex¬ 
pressed as 

( 2 ) 

Mm/Mp - <^a/Wp 

where = (318 334 539 ± 10) x 10“® is the 

ratio of the magnetic moments of the muon to 
the proton, which has been measured with a pre¬ 
cision of 30 ppb [4]. To avoid experimenter biases, 
a “blind analysis” strategy is pursued, separating 
the Wa and Wp analyses with secret offsets which 
are only revealed when both analyses are com¬ 
plete and internally consistent. 

The beam used for the experiment originates 
at the Alternating Gradient Synchrotron (AGS) 
which every 2.5 s delivers 40-60 x 10^^ protons at 
24 GeV/c onto a nickel target. Each proton spill 
is composed of 12 bunches with a width of about 
50 ns and a separation of 33 ms. Downstream of 
the target, pions at 3.1 GeV/c are selected into a 
72 m long straight beam line where about half of 
them decay into muons. Because of parity viola¬ 
tion in the pion decay, the selection of forward¬ 
going muons leads to a polarisation of about 96 %. 
Muons at the magic momentum are selected and 
injected into the storage ring through a hole in 
the yoke of the dipole magnet whose 1.45 T field 
is locally cancelled by a DC super-conducting in- 
flector magnet [5]. To move the muons onto the 
central orbit, a kick of about llmrad is given 
by a pulsed kicker magnet [ 6 ]. The continuous 
superferric ‘C’-shaped storage ring magnet [7] is 
excited by superconducting coils. The muon stor¬ 
age region has a 9 cm diameter cross-section and 
a central radius of 7.112 m, corresponding to a cy¬ 
clotron period of 149.2 ns for muons at 7 = 29.3. 

A vertical air gap between pole and yoke decou¬ 
ples the yoke and pole pieces, which are fabricated 
from high quality steel, and allows the insertion of 
iron wedges to improve the field homogeneity by 


compensating the quadrupole field components. 
The four edge shims, 5 cm wide and about 3 mm 
high, are the main tool for reducing field vari¬ 
ations over the beam cross-section. Surface coils 
glued to the pole pieces are used to further reduce 
the inhomogeneity of the field. 

The field inside the storage region is mapped 
twice a week using a hermetically sealed cable-car 
with a matrix of 17 NMR probes moving on rails 
in the vacuum beam pipe and measuring a trans¬ 
verse field map about every 5 mm. The probes 
in the trolley are calibrated in place relative to a 
standard H 2 O probe for which the calibration ra¬ 
tio u;p(H 20 )/wp(free) between the precession fre¬ 
quencies of protons bound in water and free pro¬ 
tons is known to 50 ppb [ 8 ]. 

During data taking, an array of 375 NMR 
probes embedded in the top and bottom plates of 
the vacuum chamber is used to monitor magnetic 
field variations between trolley measurements and 
to stabilise the field with a feed-back loop to the 
main magnet power supply [9]. 

Two independent analyses determined (ujp), 
averaged over the muon distribution. Their re¬ 
sults agree within 0.05 ppm. The systematic 
uncertainties for {up) are summarised in Ta¬ 
ble 1. As a final result, the value (wp)/(27r) = 
61 791 595(00)(15)Hz (0.2ppm) was obtained. 
The improvement from 0.4 ppm systematic error 
in 1999 to 0.24 ppm in 2000 comes mainly from 
the better field homogeneity (Figure 1) which 
was achieved by replacing the old inflector whose 
super-conducting fringe-field shield had a flux 
leak. 

The decay positrons from /i’*' —> range 

in energy from 0 to 3.1 GeV, and are detected 
with 24 lead/scintillating-fiber calorimeters [10] 
placed symmetrically around the inside of the 
storage ring. The arrival times and energies of 
the positrons are determined from the calorime¬ 
ter pulses whose full shapes are sampled by 
a 400MHz waveform digitiser (WED). A laser 
and light-emitting-diode (LED) system is used to 
monitor potential time and gain shifts. 

Because of parity violation in the weak muon 
decay, in the muon rest frame positrons are pref¬ 
erentially emitted along the muon spin direc¬ 
tion. Since in the lab frame forward positrons 
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Figure 1. Magnetic field map of the storage ring 
cross-section averaged over azimuth. The contour 
lines have a distance of 0.5 ppm. 

Table 1 


Systematic uncertainties for the uop analysis. 


Source of errors 

Size [ppm] 

Absolute calibration of standard probe 

0.05 

Calibration of trolley probes 

0.15 

Trolley measurements of central 


azimuthal average field 

0.10 

Interpolation with fixed probes 

0.10 

Uncertainty from muon distribution 

0.03 

Others'^ 

0.10 

Total systematic error on uip 

0.24 


^ higher multipoles, trolley temperature and volt- 
age response, eddy currents from kickers, and time- 
varying stray fields. 



are boosted to high energies, the muon spin pre¬ 
cession frequency modulates the decay positron 
count rate N{t) if a lower energy threshold is ap¬ 
plied: 

N{t) = No e~~ [I + A cos (cuat + 4>a)] , (3) 

where yr = 64.4 /is is the dilated muon lifetime. 
For an energy threshold of 2 GeV, the asymme¬ 
try A is about 0.4. Figure 2 shows the sum of the 
decay positron time spectra observed by all de¬ 
tectors within a time range of 805 ^s or roughly 
12 muon lifetimes. The total number of positrons 
recorded later than 45 /iS after injection is about 
4 X 10®. Error bars are drawn on all points but 
are only visible at very late times because of the 


huge statistics of up to 10^ entries per 149 ns time 
bin. Given this high number of events, the simple 
parametrisation in Eq. (3) proved not to be ad¬ 
equate for fitting the spectrum which is affected 
by several significant perturbations. 


4Billion Positrons with E >2 GeV 



Figure 2. Positron time spectrum from the 2000 data 
set with all detectors combined. 

We shall concentrate here on the leading sys¬ 
tematic effect caused by coherent betatron oscilla¬ 
tions - oscillations of the beam as a whole. Since 
the inflector aperture is smaller than the storage 
ring aperture, the phase space for betatron oscil¬ 
lations is not filled, which leads to a radial modu¬ 
lation of beam width and centroid. Looking from 
fixed detector positions, these oscillations have 
a frequency of 466 kHz, approximately given by 
/CBO (1 - VI - ri)fc where fc = 6.7 MHz is 
the cyclotron frequency and n = 0.136 is the field 
index adjusted by the electric quadrupole volt¬ 
age [3]. The value of n was chosen far away from 
physical resonance conditions which could lead to 
increased muon losses or spin flips. Since both 
the calorimeter acceptance and the energy dis¬ 
tribution of the detected positrons depend on the 
radial position of the muon decay, the parameters 
No, A and (/)a in Eq. (3) receive a time modula¬ 
tion, e.g. No{t) = No[l + gN(t)cos{2TTfcBot + 
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^^CBo)]) where gNit) describes the CBO decay 
due to the muon momentum spread and higher 
E and B field multipoles. The time constant is 
typically 100/rs. The modulation No{t) - with 
^^(0) « 1 % - was already observed and included 
in the analysis of the 1999 data set. Its parame¬ 
ters do not correlate strongly with tOa- The mod¬ 
ulations A{t) and </)a(t) are smaller - 0.1% and 
1 mrad at t = 0 respectively - and were only dis¬ 
covered with the higher statistics of the 2000 data 
set. Their dangerous effect on the time spectrum 
is the creation of interference terms with the fre¬ 
quency /cBO ~ Wo/(27r) which for our choice of 
n is very close to Ua- This parametric resonance 
can produce shifts in uja of up to 4 ppm if indi¬ 
vidual detector time spectra are fitted without 
including A{t) and (pa{t)- In the sum of all de¬ 
tector spectra, the shifts largely cancel thanks to 
the circular symmetry of the (g-2) ring. Remain¬ 
ing effects are accounted for by including A{t) 
into the fit and assigning a systematic error for 
the impact of (pait). Including </a(t) into the fit 
turned out to be more difficult because this term 
induces strong correlations between detector gain 
changes and uJa, entailing further shifts in uja- 

Other perturbations were treated like in ear¬ 
lier data sets [11]. Pulse pileup effects were 
removed by statistically superimposing recorded 
pulses from the data themselves and thus con¬ 
structing an artificial pileup spectrum which was 
then subtracted from the untreated spectrum. 
The effects of beam debunching were eliminated 
by randomizing the start time of each fill over one 
cyclotron period. Muon losses were taken into 
account by multiplying the function in Eq. (3) 
with an extra loss term. AGS background due 
to erroneous proton extraction during the muon 
storage period, which contributed 0.1 ppm uncer¬ 
tainty to the 1999 result, was largely eliminated 
by installing a sweeper magnet in the beam line. 

Four independent uJa analyses with different ap¬ 
proaches to take systematic effects into account 
were completed. Their results agreed within 
0.4 ppm - as compared to statistically allowed 
variations of 0.5 ppm - and were combined to 
LOa/i2Tr) = 229 074.11(14)(7)Hz (0.7ppm). This 
number contains a correction of -1-0.76(3) ppm for 
residual effects of the electric field on muons with 


7 29.3 and for deviations from Eq. (1) due 

to vertical beam oscillations (i.e. (3 ■ B ^ 0). 
The combined systematic errors listed in Table 2 
account for the correlations between the results 
from the individual analyses. 

Table 2 


Systematic uncertainties for the Wo analysis. 


Source of errors 

Size [ppm] 

Coherent betatron oscillations 

0.21 

Pileup 

0.13 

Detector gain changes 

0.13 

Lost muons 

0.10 

Binning and fitting procedure 

0.06 

Others! 

0.06 

Total systematic error on Wa 

0.31 


^ AGS background, timing shifts, E field and vertical 
oscillations, beam debunching and randomisation. 


3. Result from the Data Set of 2000 and 
Comparison with Theory 

After completion of the lOa and Wp analyses, 
was calculated according to Eq. (2). The result is 
a^+ = 11659 204(7)(5) x 10-i° (0.7ppm) [12]. 
It agrees well with older measurements (Fig¬ 
ure 3). The new experimental world average, 
a^+ = 11659 203(8) x 10“^° (0.7ppm), is domi¬ 
nated by the new result which has about half the 
uncertainty of previous measurements. 

The Standard Model prediction for [14] can 
be written as 

a^(SM) = a^(QED)-|-a^(weak)-|-ap(hadronic)(4) 

with a^(QED) = 11 658 470.57(0.29) x 10-^° and 
a;^(weak) = 15.1(0.4) x 10“^°. The hadronic con¬ 
tribution cannot be calculated from first prin¬ 
ciples at this time because it is dominated by 
low-energy interactions. The first-order hadronic 
vacuum polarisation contribution, ay^(had, 1), can 
be determined from measured e“'"e“ annihilation 
cross-sections over all energies using a dispersion 
relation. It can also be related to hadronic r 
decay. These calculations are still under theo¬ 
retical investigation. Higher-order contributions 
are given by a^(had, 2) = —10.0(0.6) x 10“^° 
and a^(had, light-by-light) = -1-8.6(3.2) xl0“^°. 
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Figure 3. Recent experimental and theoretical val¬ 
ues of The labels for the theoretical predictions 
represent the evaluation of a^(had, 1) (see text). 


For the latter, a sign error was recently cor¬ 
rected. Figure 3 shows the theoretical predic¬ 
tions for using recent evaluations of a^(had, 1). 
DH’98 [13] is the value used for the comparison in 
our PRL [12]. In this evaluation, a;_j(had, 1) uses 
data from both e“'"e“ annihilation and r decay. 
Later, new e+e” data from Novosibirsk [15] gave 
rise to a new calculation of a^(had, 1) which does 
not agree any more with the value obtained from 
T decay data (DEHZ’02: [14]). Using the e+e“- 
based result, a^(had, 1) = 684.7(7.0) x 10“^°, 
one obtains a total theory prediction a^(SM) = 
11659169.1 (7.8) x 10-^° (0.7ppm). Using the 
T-based result, o^(had, 1) = 701.9(6.1) x 10“^°, 
one obtains a total theory prediction a^(SM) = 
11659186.3(7.1) x 10-i° (0.6ppm). The devi¬ 
ations of the two evaluations from the experi¬ 
mental result correspond to 3.0 and 1.6 standard 
deviations, respectively. Hence, no unambigu¬ 
ous statement about new physics can be made 
at present. 

4. Outlook 

In the year 2001, the experiment was performed 
with negative muons and with different field in¬ 
dices n moving the CBO frequency away from the 
parametric resonance. The resulting data set of 
about 3 X 10® electrons is currently being ana¬ 


lysed. It will provide a test of CPT invariance 
and - if CPT holds - an improved combined value 
of o^±. However, in order to achieve the de¬ 
sign goal of 0.35 ppm statistical uncertainty, ad¬ 
ditional 6 X 10® events are needed. A new run is 
planned, but at present funding is uncertain. 
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